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ABSTRACT 



f 



Context. The number of known Be/X-ray binaries in the Large Magellanic Cloud is small compared to the observed population of 
the Galaxy or the Small Magellanic Cloud. The discovery of a system in outburst provides the rare opportunity to measure its X-ray 
properties in detail. IGR J054 14-6858 was serendipitously found in outburst with the Swift satellite in August 201 1. 
Aims. In order to characterise the system, we analysed the data from a follow-up XMM-Newton target of opportunity observation and 
investigate the stellar counterpart with photometry and spectroscopy. 

Methods. We modelled the X-ray spectra from the EPIC instruments on XMM-Newton and compared them with Swift archival 
data. In the X-ray and optical light curves, we searched for periodicities and variability. The optical counterpart was classified using 
spectroscopy obtained with ESO's Faint Object Spectrograph at NTT. 

Results. The X-ray spectra as seen in 201 1 are relatively hard with a photon index of ~0.3-0.4 and show only low absorption. They 
deviate significantly from earlier spectra of a probable type II outburst in 2010. The neutron star spin period of P spin = 4.4208 s was 
discovered with EPIC-pn. The /-band light curve revealed a transition from a low to a high state around MJD 54500. The optical 
counterpart is classified to B0-1 Hie and shows Ho- emission and a variable NIR excess, vanishing during the 2010 outburst. In the 
optical high state, we found a periodicity at 19.9 days, probably caused by binarity and indicating the orbital period. 

Key words, galaxies: individual: Large Magellanic Cloud - stars: emission-line, Be - stars: neutron - X-rays: binaries 



1. Introduction 

Depending on the donor star, high mass X-ray binaries 
(HMXBs) are divided into super-g iant systems and Be/X-ray bi- 
naries (BeXRBs, e.g. Reig 201 ll) . In the latter case, a not well 
understood mechanism causes matter ejection of the Be star in 
the equatorial plane, leading to the build up of an equatorial de- 
cretion disc around the Be star (Okazaki2001). This disc domi- 
nates the emission of the system in infrared and some emission 
lines like Ha. The variability of this emission points to the insta- 
bility of such discs. Due to a supernova kick, a neutron star (NS) 
can have an eccentric orbit around the Be star. During periastron 
passage the NS can accrete matter from the decretion disc, caus- 
ing a so-called type I X-ray outburst, enduring several days at 
typical luminosities of 10 36 erg s _1 . During disc instabilities, the 
NS can accrete a large fraction of the decretion disc, resulting in 
type II outbursts with luminosities up to 10 37 erg s _1 for several 
weeks. 

The transient behaviour of Be/X-ray binaries and the wide 
extent of the Large Magellanic Cloud (LMC) on the sky, which 
imposes a low observational coverage by X-ray missions, com- 
plicate the discovery and investigation of Be/X-ray binaries in 
this galaxy. In contrast t o that, the SMC was monitored with 
RXTE for about 14 years dGalache et al.ll2008h . Therefore, only 
nine HMXB X-ray pulsars are known to date in the LMC, which 



inhibits a statistical comparison of this sample with those of the 
Galaxy and the Small Magellanic Cloud (SMC). In the Galaxy 
and the SMC ~66 and ~55 HMXB pulsars are known, respec- 
tively. A major fraction o f the pulsars is found in Be/X-ray bina- 
ries (e.g. ICoe et al1l2010h . Population studies of these systems 
are important to understand the stellar evol ution, as th ey e.g. 
allow to estimate sup ernova kick velocities (|Co e 2005) or th e 
star formation history (lAntoniou et alfe OlO: Min eo et alJl201 ll) . 
Recently, a bimodal NS spin period distribution for the Galactic 
and SMC samp les has been assoc iated with two different types 
of supernovae dKnigge et al.l l20l"lh . To enable such statistical 
studies for the LMC, it is necessary to successively build up a 
larger sample of X-ray measurements of pulsars in outburst. 

In 2010, IGRJ05414-6858 was discovered serendip- 
itously within INTEG R AL observations of SN 1987 A 
(Grebenev & Lutovinov I20T0I) and later localised 
(Lutovinov & Grebenev 120101) and identified as Be-X-ray 



* Based on observations with XMM-Newton, an ESA Science 
Mission with instruments and contributions directly funded by ESA 
Member states and the USA (NASA) 



binary dRauetal.1 1201 Oh with Swift and GROND follow-up 
observations. In 201 1, Swift performed a UV survey of the LMC 
(PI: S. Immler). This provided a shallow coverage of this galaxy 
with the Swift X-ray telescope (XRT) and allowed the detection 
of bright X-ray transients. In an observation on 2011 Aug. 
5 an o utburst of IGR J05414-6858 was detected dSturm et all 
1201 lbl) . which allowed us to request an XMM-Newton target of 
opportunity (ToO) observation. 

In this study, we report our analysis of the XMM-Newton 
observation of IGR J05414-6858. The detection of the NS spin 
period adds the tenth X-ray pulsar in the LMC sample and a de- 



1 



Sturm et al.: The Be/X-ray binary IGR J05414-6858. 
Table 1. X-ray observations of IGRJ05414-6858. 



Observation 


ObsID 


Date 


Time 


Instrument 


Mode" 


Offax fc 


Net Exp 


Net Cts'' 


R sc 








(UT) 






['] 


[ks] 




["] 


XMM 2001 


f\f\f\ A A 1 A 1 A 1 

0094410101 


2001-10-19 


11:49 - 20:30 


EPIC-pn 


ff-medium 


8.2 


8.2 


<26 


30 








17:09 - 20:29 


EPIC-MOS1 


ff-medium 


7.2 


11.7 


<5 


30 








17:09 - 20:29 


EPIC-MOS2 


ff-medium 


6.9 


11.7 


<14 


30 


XMM 2011 


0679380101 


2011-08-13 


07:45 - 14:07 


EPIC-pn 


ff-thin 


1.1 


7.2 


2018 


21 








07:22 - 14:07 


EPIC-MOS 1 


ff-medium 


1.2 


9.4 


863 


29 








07:22 - 14:07 


EPIC-MOS2 


ff-medium 


1.2 


9.4 


954 


31 


Swift 2010 a 


00031745001 


2010-06-25 


05:35 -09:10 


XRT 


pc 


5.3 


3.9 


221 


40 


Swift 2010 b 


00031745002 


2010-06-30 


01:14-23:59 


XRT 


pc 


2.0 


5.2 


497 


50 


Swift 2011 a 


00045428001 


2011-08-05 


15:28-22:26 


XRT 


pc 


10.5 


3.3 


40 


40 


Swift 2011 b 


00031745003 


2011-08-09 


23:53-04:55 


XRT 


pc 


1.8 


3.7 


120 


40 


Swift 2011 c 


00031745004 


2011-08-12 


00:03-22:53 


XRT 


pc 


2.4 


4.1 


100 


40 


Swift 2011 d 


00031745005 


2011-08-20 


00:25 - 16:53 


XRT 


pc 


3.1 


2.1 


<6 


40 


Swift 2011 e 


00031745006 


2011-08-24 


02:44 - 04:29 


XRT 


pc 


2.5 


1.2 


<5 


40 



Notes. <D) Observation setup: full-frame mode (ff) and photon counting mode (pc). For XMM-Newton, also the filter is given. {b) Off-axis angle 
under which the source was observed. (r) Net counts as used for spectral analysis in the (0.2 - 10.0) keV band for XMM-Newton and in the (0.3 - 
6.0) keV band for Swift. {d) Radius of the circular source extraction region. 



tailed X-ray spectral and temporal analysis allows a characteri- 
sation of the system. We compare our new results to those from 
archival Swift data and discuss complementary optical data to 
characterise the optical counterpart and the circum-stellar disc. 



2. Observations and data reduction 

2.1. XMM-Newton 

The XMM-Newton dJansen et al.ll200ll) ToO observation was 
performed on 2011 Aug. 13. The source was obs erved on-axis, 
placed on CCD4 of EPIC-pn dStriider et al.ll200ll) and CCD1 of 
both EPIC-MOS (Turne r et al Il200ll) detectors. We used XMM- 
Newton SAS ll.O.cQto process the data. Unfortunately, the ob- 
servation was affected by an increased background caused by 
soft protons. During the first -11.5 ks of the ~17 ks observa- 
tion, the background was at a moderately elevated level, allow- 
ing the selection of time intervals where the background rate in 
the (7.0-15.0) keV band was below 50 cts ks -1 arcmin -2 for 
EPIC-pn and below 4 cts ks 1 arcmin" 2 for EPIC-MOS. The de- 
tailed observation setup is recorded in TableQ] Here, we also list 
an XMM-Newton observation from 2001 covering the position 
of IGR J05414-6858, in which the source was not detected. We 
used this observation to derive an upper limit for the flux. For 
EPIC-pn, we used single and double pixel events and single to 
quadruple events in the case of EPIC-MOS, all having FLAG=Q. 
Background events were selected from a point source free area 
on the same CCDs as the source. Source events were extracted 
from a circle, with radius optimised for the signal-to-noise ratio 
by the SAS task eregionanalyse. We created spectra and re- 
sponse matrices with especget and used a binning to have at 
least a signal to noise ratio of 5 for each bin. For time series, 
the photon arrival times were randomised within the CCD frame 
time and calculated for the solar system barycentre. 



1 Science Analysis Software (SAS), http://xmm.vilspa.esa.es/sas/ 



2.2. Swift 

We re-analysed archival Swift/XKT observations. The spectra 
were created by using the ftoo]Q xselect to select events in 
the cleaned level 3 event files within a circle, placed on the 
source with radii given in Table [TJ Background spectra were 
created from a circular extraction region with radius of 200". 
The spectra were binned to have >20 cts bin -1 . The ancillary re- 
sponse files were calculated with xrtmkarf . Swift observations 
of IGR J05414-6858 including non-detections, are also listed in 
Table [TJ 



2.3. INTEGRAL 

For INTEGRAL/ISGRI we used the data taken in May and 
Jun 2010, i n which the source was disco vered in the (20^40) 
keV band dGrebenev & Lutovinovl |20T5T) . We selected point- 
ings within 4.8 degrees around the source, so that the source 
was in the fully coded field of view of ISGRI and excluded 
pointings with timing problems (non-zero TIMECORR). We 
had 135 pointings and 470 ks exposure after the selection. 
For the data analysis we used the official INTEGRAL soft- 
ware "Offline Science Analysis" (OSA) version 9.0 provided by 
the INTEGRAL Science Data Center (ISDCfl Because there 
is no straight forward way of obtaining light curves on short 
time scales for weak sources in the standard pipeline, we did 
the timing analysis without deconvolution. We ran the standard 
analysis from the COR to the DEAD level, created the Photon 
Illumination Function map for IGRJ05414-6858 with the stand- 
alone OSA tool ii_pi£, then extracted single events from good 
pixels (SELECT_FLAG==0) in the energy range of 15 keV to 
300 keV with evts_extract, applying the barycentric correc- 
tion during the course of events extraction, then selected the 
typeO events with PIF==1 with f tools. 



2 http://heasarc.nasa.gov/ftools/ 

3 http://www.isdc.unige.ch/integral/analysis 
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2.4. NIR, optical, and UV photometry 

Optical photometry of IGR J05414-6858 was obtained with 
OGLE, GROND and Swift/UWOT. The optical counterpart was 
monitored regularly during the Optical Gravitationa l Lensing 
Experiment (OGLE) of phase III dUdalski et al.ll2008l) between 
October 2001 and April 2009 in the /-band. The source identifi- 
cation is OGLEIIILMC175.4.21714. 

The Gamma-ra y Burst Optical Near-ir Detector (GROND; 
iGreiner et af]|2008h at the MPG/ESO 2.2m telescope in La Silla, 
Chile, observed the source at three epochs in June 2010 and 
January 2012. Preli minary resu l ts of th e 2010 observations were 
already presented in iRau etafl d2010l) . GROND is a 7-channel 
imager that observes in four optical and three near-IR chan- 
nels simultaneously. The IGR J05414-6858 data were reduced 
and analysed with t he standard tools and methods described in 
iKriihler et al.l d2008l) . The photometry was obtained using point- 
spread-function (PSF) fitting taking into account the contami- 
nation from the two nearby sources (see Fig. Q]). Calibration 
was performed against observations of an SPS S standard star 
field (g'r'i'z!) or against selected 2MASS stars (Skrutskie et al. 
120061) (JHK S ). This resulted in lcr accuracies of 0.04 mag (g'z'X 
0.03 mag (r'/'), 0.05 mag (JH), and 0.07 mag (K s ). 

The Ultraviolet/Optical Telescope (UVOT) aboard Swift has 
three optical (v b u) and three UV filters (uvwl uvm2 uvw2). 
IGRJ05414-6858 was observed in all UVOT filters during the 
pointed X-ray observations. For the LMC UV-survey observa- 
tion 00045428001, the source is not in the field of view of the 
UVOT. Photometry was carried out on pipeline processed sky 
images downloaded from t he Swift data cen trjj following the 
standard UVOT procedure dPoole et al.ll2008h . 

2.5. Optical spectroscopy 

Optical spectroscopy was taken with the ESO Faint Object 
Spectrograph (EFOSC2) mounted at the Nasmyth B focus of 
the 3.6m New Technology Telescope (NTT), La Silla, Chile on 
the nights of 2011 December 8 and 10. The EFOSC2 detector 
(CCD#40) is a Loral/Lesser, Thinned, AR coated, UV flooded, 
MPP chip with 2048x2048 pixels corresponding to 4.1'x4.1'on 
the sky. The instrument was in longslit mode with a slit width 
of 1.5". Grisms 14 and 20 were used for blue and red end 
spectroscopy respectively. Grism 14 has a wavelength range of 
/L13095-5085 A and a grating of 600 lines mirr 1 and a disper- 
sion of 1 A pixeL 1 . The resulting spectra have a spectral reso- 
lution of ~ 5 A. Grism 20 is one of the two new Volume-Phase 
Holographic grisms recently added to EFOSC2. It has a smaller 
wavelength range, from 6047-7147 A, but a superior disper- 
sion of 0.55 A pixeL 1 and 1070 lines mrrT 1 . This produced 
an almost identical spectral resolution for our red end spectra, 
~ 10 A. Filter OG530 was used to block second order effects. The 
data were reduced using the standard packages available in the 
Image Reduction and Analysis Facility I RAF. Wavelength cal- 
ibration was implemented using comparison spectra of Helium 
and Argon lamps taken throughout the observing run with the 
same instrument configuration. The spectra were normalized to 
remove the continuum and and a redshift correction applied cor- 
responding to the recession velocity of the LMC (-280 km s _1 , 
iPaturel et alJl2002h . 



4 http://w w w. swift, ac.uk/swift_portal 
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RA (J2000) 

Fig. 1. GROND r'-band finding chart. Lines mark the counter- 
part of IGR J05414-6858. In the zoom-in, the XMM-Newton po- 
sition is marked by a white circle with radius of the lcr position 
uncertainty of 0.52". The astrometric solution accuracy of the 
GROND image is 0.47" in RA and 0.21" in Dec. 

3. Analyses and results of X-ray data 

3. 1 . X-ray coordinates 

We created X-ray images from all three EPIC cameras in 
the XMM-Newton standard energy sub-bands. A simultane- 
ous source detection was performed on these images with 
edetect.chain. The best-fit source position is RA (J2000) = 
05 h 41 m 26 s .62 and Dec (J2000) = -69°01'23'/0. The lcr uncer- 
tainty of the position is 0.52", where we assume a systematic 
error of 0.5", which is quadratically added to the statistical er- 
ror. The angular separation to the optical counterpart is 0.52" 
for the 2 MASS position and 0.69" for the GROND position 
("st ar A. IRau et al.l l2010h . The distance to the Swift position 
of iLutovinov & Grebenevl (l2010b is 2.0" with an uncertainty in 
the Swift measurement of ~3". A finding chart obtained from 
GROND data is presented in Fig. Q] The white circle in the 
zoom-in gives the XMM-Newton position. The improved X-ray 
coordinates from the XMM-Newton observation further confirm 
the identification of the X-ray source with the optical counter- 
part. 

3.2. Spectral analysis 

Spectral analysis was performed with xspec (lArnau d1 fl996h 
version 12.7.0. The three XMM-Newton EPIC spectra were 
fitted simultaneously and we included constant factors in the 
models to consider instrumental differences. For all models, 
we obtained consistent values of Cmosi = 1.11 ± 0.08 and 
C M os2 = L10 + 0.08 relative to EPIC-pn (C pn = 1). Therefore, 
the fluxes for all instruments are consistent within uncertain- 
ties, as EPIC-MOS is known to derive ~5% higher values com- 
pared to EPIC-pn. All other model parameters were forced to 
be the same for all instruments. The spectra are well described 
by an absorbed power-law. The photoelectric absorption was 
modelled by a fixed G alactic foreground column density of 
A^H,Gai = 6 x 10 20 c m" 2 dDickev & Lockmanll 19901) with abun- 
dances according to I Wilms et al. (2000|). An additional column 
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density A^h.lmc was determined in the fit. It accounts for the in- 
terstellar medium of the LMC and source intrinsic absorption 
and the abundances were set to 0.5 for elements heavier than 
helium dRussell & Dopitall 19921) . 

We further tested the spectra for the existenc e of typical fea- 
tures of BeX RBs. A possi ble soft excess (e.g. Eg er & Haberll 
2008; iHickox et alj 120041) was modelled by a black-body 
(Model: PL+BB) or a multi-temperature disc black-body model 
(PL+DiscBB). The additional model component only improves 
the fit marginally, but demonstrates model-dependent uncertain- 
ties for the power-law parameters. The f-test probability for these 
additional components is around 13%. We use the black-body 
model to derive an upper limit for a soft component contribution. 
A soft excess can account only for <1% of the detected flux and 
up to 22.0% of the unabsorbed luminosity in the (0.2-10.0) keV 
band. The derived radius of the emission region from the black 
body is too large for an NS. With the disc model we obtain a sim- 
ilar inner radius f or an inclination of = (/?;„ oc (cos ©) ^ 2 ). 
Also according to Hick ox et alj (12004), we obtain an inner disc 
radius of R[ D = (Lx/4-jtctT 4 ) 1 ^ 2 = 52 km. Fluorescent iron line 
emission at 6.4 keV was modelled with a Gaussian line, having 
a fixed central energy and no measurable broadening. Also here, 
we receive only a marginal improvement of the fit and a line flux 
of 2.8 + 2.6 x 10~ 6 photons crrT 2 s . The spectrum and best-fit 
model are presented in Fig|2]and the best-fit results are listed in 
Table [2] All uncertainties and limits correspond to a 90% confi- 
dence level (Ax 2 = 2.71). 

To investigate the spectral variability, we compared the EPIC 
spectra to the archival Swift spectra. The lower statistics only al- 
lowed to fit the power-law model. For the two spectra in 2010, 
jVh.lmc and T are consistent within errors. Also, no significant 
evolution is found in the three Swift spectra of 201 1 . We assumed 
that the spectral shape is time-independent during the individual 
outbursts, and fitted the model simultaneously to the two 2010 
and three 2011 Swift spectra (see Fig. [3J. The fluxes were de- 
termined for all observations individually. The results are also 
listed in Table|2] There is significant spectral variability observed 
between the outbursts of IGR J05414-6858 in 2010 and 2011, 
where the recent outburst exhibits a harder X-ray spectrum, con- 
sistent with the XMM-Newton 2011 observation. If we fit the 
Swift 2010 data with the spectral model derived from the XMM- 
Newton 2011 observation (allowing a re-normalisation), the fit 
degrades from^ 2 ed = 1.03 to 1.75. Since the lower photon index 
in 2011 is seen with XMM-Newton and Swift, this is unlikely 
caused by the high background during the XMM-Newton obser- 
vation or by instrumental differences. 

3.3. Pulsations 

A strong signal at co = 0.2262 Hz and its first harmonic ap- 
peared in a Fast Fourier Transformation (FFT) of the EPIC- 
pn time series in the (0.2-10.0) keV band. The power density 
spectrum is plotted in Fig. |4] The signal is also clearly present 
in the (0.2-2.0) keV and (2.0-10.0) keV sub-bands. The pe- 
riod is not resolved by EPIC-MOS (2.6 s frame time), since the 
period is shorter than twice the frame time of the instrument, 
i.e. cj is above the Nyquist frequency. The same holds for the 
Swi ft data (2.5 s frame time ). A x 2 test, a Bayesian odds ra- 
tio dGreeory & Loredolll996l). and a Rayleigh Z 2 test for one 
harmonic (Haberl &Zavlin 2002; Buccher i et al.llT983b around 
the periodicity signal are shown in Fig. [5] All tests indepen - 
dently confirm the pulse period. Following Haberl et al. (2008), 
we used the Bayesian detection method to determine the pulse 
period and a lcr uncertainty of 4.420866(2) s on 201 1-08-13. 




0.5 1 2 5 10 

Channel Energy (keV) 



Fig. 2. EPIC-pn (black), EPIC-MOS 1 (red), EPIC-MOS2 
(green) spectra, together with the best-fit po+bb+Fe model 
(solid line) and its individual components: power-law (dashed), 
black-body (dashed-dotted) and Fe line (dotted). The lower 
panel shows the residuals. 



Fig. [6] shows the folded background-subtracted light curves 
from EPIC-pn in the total (0.2 - 10.0) keV band and the standard 
sub-bands (0.2 - 0.5) keV, (0.5 - 1.0) keV, (1.0 - 2.0) keV, (2.0 
- 4.5) keV, and (4.5 - 10.0) keV, where we merged the first two 
bands, to increase the statistics. Also hardness-ratio (HR) varia- 
tions are presented. HRs are defined by HR, = (Rj+i ~ R/)/(Ri+i 
+ R,) with R,- denoting the background-subtracted count rate in 
the standard energy band i (with i from 1 to 4). In the light 
curves, two narrow peaks are seen within one period having only 
small variations in energy. By modelling the (0.2 - 10.0) keV 
curve with a non-pulsating contribution and one Gaussian for 
each peak, we estimate a pulsed fraction of (48+7)% in the total 
flux, and a flux ratio for both peaks of 2.3±0.5. 

The pulsation period does not show up in the power density 
spectrum of the INTEGRAL events. Also in the 4.42 s folded 
light curve, we did not find a significant peak. This might be 
caused by orbital modulations, as the INTEGRAL observations 
cover a long time. 

3.4. X-ray flux variability 

In addition to the fluxes as measured in Sec. 13.21 we calculated 
upper limits for non-detections of IGR J05414-6858. The field 
was observed by XMM-Newton in 2001, but no source was de- 
tected by edetect_chain (analogous to Sec. 13. 11 1. Spectra were 
extracted in the same manner as described in Sec. [2] from a 
30" source region and a 50" background region. Using C statis- 
tics and the spectral shape as determined with XMM-Newton in 
2011, we derived a 90% upper limit for the flux of 9.0xl0~ 14 
erg cirT 2 s _1 . This lowest upper limit results in a variability of at 
least a factor of 100, compared to the maximum flux measured 
in the Swift 2010b observation. Prior to that, we could not find 
any corresponding X-ray detection. There is no ROSAT source 
within 1' listed in the literature. 

The Swift monitoring of the recent outburst, determined the 
turn-off between MJD 55786 and 55793. For the two non- 
detections, we extracted spectra in the same way as above and 
determined 90% confidence upper limits of 9.5 x 10~ 13 erg cm" 2 
s _1 and 2.3 x 10~ 12 erg crrT 2 s _1 . The long term evolution of the 
system is presented in Fig. [7] 
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Table 2. Spectral fit results. 



Observation 


Model" 


A'h.lmc 
[10 21 cm- 2 ] 


r 


[eV] 


[km] 


EW Fe 

[eV] 


Flux c 
[10~ 12 erg cm -2 s~'] 


[10 36 ergs" 1 ] 


^red 


dof 


XMM 2011 


PL 

PL+BB 

PL+BB+Fe 

PL+DiscBB 


<0.89 
4.98£« 
4.98+J 6 , 
3.681" 


32 

"•-'^-0.05 

40 +0 07 

041 +0.07 


112+g 

lu+ 49 

128 + ff 

-33 


49 + ! 34 

-44 


78$ 


3 01 

2 98 +0 20 

9 Q7+0.37 
z ""-047 

2.98 + °i? 

-0.33 


0.91 
0.98 
0.98 
1.01 


0.92 
0.91 
0.88 
0.91 


118 
116 
115 
116 


Swift 2010 a 
Swift 2010 b 


PL 


1 n+2.0 
1 - y -l.4 


0.78^ 








5 6 +u 
9 2 +1 3 


1.77 
2.91 


1.03 


31 


Swift 2011 a 
Swift 201 1 b 
Swift 2011 c 


PL 


<2.6 


0.25333 








5 4 +26 

5 1 +U 

3 5 +98 

-3.1 


1.62 
1.51 
1.04 


1.18 


8 



Notes. <n> For definition of spectral models see text. <6) Radius of the emitting area (for BB) or inner disc radius for an inclination of = (for 
DiscBB, for the definition see text). (c) Observed flux in the (0.2-10.0) keV band, derived by integrating the best-fit model. Source intrinsic 
X-ray luminosity in the (0.2-10.0) keV band corrected for absorption and assuming a distance of the source of 50 kpc. 
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Fig. 3. Swift spectra of IGR J05414-6858 from 2010 (top) and 
201 1 (bottom) with best-fit power-law model. Lower panels give 
the residuals. 



The variability during the XMM-Newton observation in 
201 1 was at a moderate level. We created a background corrected 
light curve, merged using all EPIC instruments and binned to 70 




mm 



J'k'Hin 1 ■| J l l lryill Indian 

0.01 0.1 1 

Frequency (Hz) 

Fig. 4. Power density spectrum of IGR J05414-6858 for the 
EPIC-pn time series in the (0.2-10.0) keV band. The best-fit fre- 
quency of a) = 0.2262 Hz and its first and second harmonics are 
marked with dashed lines. 



s, corresponding to 30 cts bin 1 on average. Ax 2 test of this light 
curve against a constant resulted in^- 2 /dof = 130/121. 



4. Analysis and results of optical data 

Using GROND, the counterpart of the S wift X-ray source was 
resolved into three point sources. Star A of lRau et all (2010) was 
originally classified as of spectral type Bl-2 III from GROND 
and Swift/WOT photometry and is the most likely counterpart 
within the improved XMM-Newton X-ray position uncertainty 
ofIGRJ05414-6858. 



4.1. Photometry 

The OGLE III /-band light curve of this star is presented in 
Fig. [7] It exhibits two different brightness states with a transi- 
tion phase in between. Before MJD 54400 (left dotted line in 
Fig- Eh the /-band emission is at (15.51-15.43) mag, followed 
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Fig. 7. /-band light curve of OGLEIIILMC175.4.21714 (upper 
panel) compared to the X-ray fluxes with upper limits marked 
by arrows (lower panel). Dotted lines separate the optical low, 
transition, and high state. 



Fig. 5. Top: x 2 test for persistence of the EPIC-pn light curve, 
around trial frequencies between 0.2252 and 0.2272 Hz. Middle: 
Frequency dependence of the Bayesian odds ratio. Bottom: 
Rayleigh Z 2 test. 
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Fig. 6. Left: X-ray pulse profile of IGR J05414-6858 in various 
energy bands from the EPIC-pn time series. The pulse profiles 
are background subtracted and normalised to the average net 
count rate of 3.0, 4.6, 6.8, 8.3 and 22.4 xl0~ 2 cts s -1 from top 
to bottom. Right: Hardness ratios as a function of pulse phase 
derived from the pulse profiles in two neighbouring standard en- 
ergy bands. 



by an impetuous increase to the high state. After MJD 54650 
(right dotted line in Fig. 13, the /-band emission is in between 

(15.24- 15.13) mag. 

A Lomb-Scargle dLomhlll97d LScarglelll982l) periodo gram 
of the low and high state of the OGLE III light curve is shown in 
Fig. [8] The high state reveals a periodicity with 19.9 days. The 
folded light curve of the high state is presented in Fig. [9] Dotted 
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Fig. 8. Lomb-Scargle periodogram of the /-band of 
OGLEIIILMC175.4.21714 for the low state (upper panel) 
and high state (lower panel). 



and dashed lines mark the phases of X-ray detections and non 
detections of 201 1, respectively. 

The Swift/U VOT photometry is presented in Table [3] Within 
uncertainties most values are constant. For the 2010 measure- 
ments, we see an indication for a flux decrease in the uvm2 and 
uvw2 magnitudes of -0.13+0.04 mag and -0.17+0.04 mag, re- 
spectively. Other magnitude in 201 1 are constant, except a pos- 
sible short increase in the UV on 201 1-08-13 between 22:49 and 
23: 12 by ~0. 1 mag. These are of the same order as the variations 
observed by OGLE during high and low state. 

By comparing the averaged magnitudes of 2010-06-25 and 
those of 201 1, we find a stronger flux increase as observed with 
OGLE III in 2008, with Av = -0.36 mag, Ab = -0.30 mag, 
Aw = -0.31 mag, Awvwl = -0.20 mag, ts.uvm.2 = -0.20 mag, 
and Amvw2 = -0. 1 8 mag. This indicates a further transition from 
low to high state between June 2010 and August 201 1. 

The GROND magnitudes are summarised in Table |4] Due to 
pointing constraints, the first observation was performed only in 
the J, H, and K s bands, while the remaining epochs cover all 
seven channels. Note the increase in magnitude between 2010- 
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Fig. 9. /-band light curve from OGLE III in the high state 
convolved with 19.898 days. Phase = corresponds to MJD 
54640. Dotted lines give the phase of the beginning of the X- 
ray observations in 201 1 . Dashed lines mark the Swift/XRT non- 
detections. 
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Fig. 12. Spectrum of IGR J045 14-6858 in the wavelength range 
/L16000-7000A with the NTT on 201 1-12-10. The spectrum has 
been smoothed with a boxcar average of 3, normalized to remove 
the continuum and shifted by -280 km s . 
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Fig. 10. UV-near-IR SED composed of Swift/UVOT observa- 
tions obtained on 2010 June 25th and GROND data taken on 
2010 June 26th. The solid line shows the simplified best fit 
thermal model indicating that the UVOT photometry is consis- 
tent with the B0-1 stellar classification derived from the spec- 
troscopy. 



06-26 and 2010-06-29 of AJ = (0.25 ± 0.07) mag, AH = 
(0.45 + 0.07) mag, and AK S = (0.42 + 0.11) mag within three 
days. This affirms the observed drop of the flux as observed with 
Swift/UVOT. 

The fit to the spectral energy distribution composed from the 
first GROND epoch together with averaged UVOT photometry 
from 2010 June 25th, is shown in Fig.[l0] Here, all magnitudes 
were corrected for G alactic re ddening of Eg y = 0.075 mag 
dSchlegel et al.ll 19981) using the lCardelli et all d 19891) extinction 
law and for the LMC-intrin sic red dening of Eb v = 0.15 mag 
(see Sec. 14.21 1 using the iPeil ([1992) extinction law. The data are 
best fit with a hot (« 31 000 K) thermal spectrum, consistent 
with the B0-1 III stellar classification suggested by the optical 
spectroscopy (Sec. l4.2t . We note the existence of a clear excess 
in the near-IR bands. 



4.2. Spectral classification 

OB stars in our own galaxy are classified using the ratio of 
certain metal and helium lines (Walborn & Fitzpatrick 1990) 
based on the Morgan-Keenan (MK; iMorgan et all 119431) sys- 
tem. However, this is unsuitable in lower metallicity environ- 
ments as the metal lines are either much weaker or not present. 
As such, the optical spectrum of IGR J05414-6858 was classified 
using the method developed by iLennonl d 1 9971) for B-type stars 
in the SMC and impleme nted for the SMC, LMC and Galaxy by 
dEvans et al.ll2004l2007h . This system is normalized to the MK 
system such that stars in both systems show the same trends in 
their line strengths. The lumi nosity classification method from 
Walb orn & Fitzpatrickl (l990) was assumed in this work. 

Fig. [TT] shows the unsmoothed optical spectrum of 
IGR J045 14-6858. The spectrum is dominated by the hydrogen 
Balmer series and neutral helium lines. The He i line at /14143 A 
is stronger than the He n /14200A, which means the star is later 
than type 09. Hen A4686 A is also clearly present implying 
that the optical counterpart of IGR J045 14-6858 is earlier than 
type B1.5, although the Hen /14541 A line is not visible above 
the noise level of the data. There is also evidence for the Siiv 
/14116 A line - consistent with a spect ral classification of Bl. 
However Wal born & Fitzpatrickl £1990) present spectra of B0 
type stars with clear Si lines indicating that a spectral classifica- 
tion of B0 is not ruled out by their presence. We note that there 
does not appear to be any evidence for the Siiv /14088 A line. 
Be stars are characterised by their rapid rotation velocity and 
it could be that this line is being concealed by the rotationally 
broadened HS line in close proximity. 

The luminosity class of the system was determined using the 
ratios of S iv A41 16/Hei 44121, He 1 44121/Hei i4143 and Hen 
/14686/Hei /14713. The first two ratios strengthen with decreas- 
ing luminosity class (i.e. with increasing luminosity) where as 
the latter ratio decreases with increasing luminosity. The relative 
strengths of these lines suggest a luminosity class III, making 
our spectral classification of B0-1 III consistent with that ob- 
tained photometrically with the GROND and Swift/UVOT data. 
To check the spectral identification we can compare the observed 
optical magnitudes with that predicted for a B0-1 III star in the 
LMC. Taking the faintest, least disc-contaminated V-band mag- 
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Table 4. GROND photometry. 
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AB Mag 
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H 


K s 




2010-06-26 10:51 
















15.93 ± 


0.05 


16.07 + 0.05 


16.29 ± 


0.07 


2010-06-29 10:29 


15.34 ± 


0.04 


15.60 + 


0.03 


15.82 + 0.03 


15.98 


t0.04 


16.17 ± 


0.05 


16.53 ±0.05 


16.71 ± 


0.08 


2012-01-25 04:14 


15.18 ± 


0.04 


15.41 ± 


0.03 


15.57 ±0.03 


15.72 


±0.04 


15.96 ± 


0.05 


16.28 ±0.05 


16.54 ± 


0.08 



Notes. {a) Not corrected for Galactic foreground reddening or extinction in the LMC. 
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Fig. 11. Spectrum of IGR J04514-6858 in the wavelength range /L13900-5000A with the NTT on 201 1-12-08. The spectrum has 
been normalized to remove the continuum and redshift corrected by -280 km s~' . Atomic transitions relevant to spectral classification 
have been marked. 



nitude from Table|3] of V = (15. 5 1 +0.06) mag, a distance modu- 
lus o f 18.5+0.1 dKoerwerl2009j) and a reddening of E B -v = 0.15 
mag ( ISchwering & Israe l 19911) reveals an intrinsic value for the 
star of V = (-3.5 + 0.2 ) mag. Such a value would be consis- 
tent with a B0. 5 III star (Wegner 2006) and hence confirms the 
classification deduced from the spectra reported here. We note 
that we cannot use the derived X-ray absorption reported here to 
refine the column to the star because of the large uncertainties in 
that value. 

Fig. [T2l shows the red end of the spectrum of IGR J04514- 
6858 taken near-simultaneously. The Ha equivalent width, con- 
sidered an indicator for circumstellar disc size, is relatively small 
at -3A which is consistent with the lack of H/3 in emission in 
Fig. QT| The double peaked, asymmetric nature of the line pro- 
file shows a V/R pattern consistent with global one armed oscil- 



lations (GOAO) and suggests that the circumstellar disc of the 
star is inclined to the line of sight. This is not uncommon in the 
circumstellar discs of Be stars. 



5. Discussion and conclusions 

We performed an XMM-Newton ToO observation of 
IGR J05414-6858 in August 2011, allowing us to measure 
the X-ray spectrum and discover the spin period. This adds the 
tenth known HMXB pulsar to the LMC sample and confirms 
the neutron star nature of the compact object. 

The first detected outburst in 2010 lasted for ~1 month, 
which suggests a type II outburst. The luminosity and observed 
duration of the 201 1 outburst are in agreement with a type I out- 
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burst, but since the time of the beginning of the outburst is un- 
known, we cannot exclude a type II outburst. 

Furthermore, the spectrum of the current outburst was found 
to be significantly harder than in 2010. The power-law photon 
index of the 2011 outburst of T — 0.3-0.4 is also relatively 
low, compared to the distribution known from the SMC sam- 
ple (T ~ 1). We note that also for a few BeXRBs in the SMC 
hard spectra were observed, as fo r the pulsars XTE J0103-728 
<T = 0.35 - 0.54, P spin = 6.85 s, lHaberl & Pietschl[2008h and 
XMMUJ0048 14.0- 732204 <T = 0.53 - 0.66, P sp]n = 11.87 s, 
ISturm et all 1201 lah . The first one of t hese also was detecte d 
in a type II outburst with INTEGRAL (iTownsend et al.ll2010l) . 
Moreover, for both SMC pulsars an indication of a soft excess 
with comparable emission radius was found and suggested to 
originate from the accretion disc. The system intrinsic absorp- 
tion strongly depends on the modelling (0-5 xlO 21 cur 2 ) and the 
location of the system in the LMC. The total LMC column den- 
sity along the lin e of sight at the p osition of IGR J05414-6858 is 
3.6 x 10 21 cirr 2 (iKim et al.ll2003l) . 

According to the Corbet relat ion dCorbetl Q984; 
lLavcock etall 120051: ICorbetetail 120091) . we expect the or- 
bital period of the system in the range of (1 - 100) days for 
the measured spin period of IGR J05414-6858. The periodic 
variations seen in the /-band are therefore likely caused by 
binarity. Assuming masses of 21.5M© and 1 M Q for the Be star 
and the NS, respectively, the orbital period implies a semi-major 
axis of the binary system of 0.41 AU, corresponding to ~6 
stellar radii (IVacca et al.lll996|) . The X-ray detections in 2011 
occurred during the bright phase of the folded /-band light 
curve (see dotted lines in Fig. [9j». In 2011, the X-ray luminosity 
follows the /-band emission and the non-detections were during 
low /-band emission. In the case of the 2010 outburst, the two 
Swift X-ray detections are during low /-band emission at phase 
0.80 and 0.04. This is either due to the fact, that type II outbursts 
are not correlated to the orbital phase, or that the optical period 
is not caused by binarity. 

Long-term optical va riation with different variabi lity patterns 
are typical for BeXRBs (Rai oelimanana et"ai1l201 lh . An expla- 
nation for the transition from low to high state by A/ ~ 0.3 mag 
might be the build up of a decretion disc around the Be star. The 
XMM-Newton upper limit (X-ray faint state) in 2001 was dur- 
ing the optical low state, where probably no decretion d isc was 
present. From other BeXRBs (c.f Fig. 10 of Reig 20TTI). a cor- 
relation between NIR and optical magnitudes is seen. Therefore 
it is likely that the long term variability seen in the OGLE /- 
band extends from NIR to UV. As indicated by the GROND and 
Svw/f/UVOT observations, the system still undergoes strong vari- 
ations in the NIR, optical and UV. A rapid strong drop in the 
NIR emission was observed between June 2010 26th and 29th. 
At this time, the source was still in a presumable type II X-ray 
outburst, which can be followed by a disc loss phase. Also, the 
/Tj-band magnitude was high, compared to the GROND January 
2012 observation, while J was at the same level again pointing 
to the presence of a circumstellar disc. This is further supported 
by the Ha line emission in December 2011. Unfortunately, the 
end of the 2010 X-ray outburst is not constrained. During the 
201 1 outburst, the Swift/UVOT observations do not indicate any 
strong changes in the optical. Forthcoming OGLE IV data will 
allow to extend the light curve and to confirm the periodicity. 

IGRJ05414-6858 is the tenth known HMXB pulsar in the 
LMC with Pspin = 4.4208 s. The optical counterpart was classi- 
fied as of spectral type B0-1 Hie and shows double-peaked Ho- 
emission and a variable NIR excess. A likely orbital period of 
^oib = 19.9 d was found. The two observed X-ray outbursts 



demonstrate the importance of optical monitoring during out- 
bursts, to better understand the accretion process in these sys- 
tems. To further increase the sample of explored HMXB in the 
LMC, further X-ray observations triggered during an outburst 
are necessary. 
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03:46 


15.05 


±0.03 










05:17 


15.05 


±0.05 




2011- 


-08- 


13 


22:52 


15.10 


±0.03 




2011- 


-08- 


-20 


17:11 


15.12 


±0.03 


uvwl 


2010- 


-06- 


■25 


05:58 


15.45 


±0.03 










07:41 


15.44 


±0.03 










09:24 


15.44 


±0.03 




2011- 


-08- 


10 


00:17 


15.20 


±0.03 










02:02 


15.25 


±0.03 










03:45 


15.24 


±0.03 










05:15 


15.26 


±0.03 




2011- 


-08- 


-13 


22:49 


15.31 


±0.03 




2011- 


-08- 


■20 


17:08 


15.22 


±0.03 


uvml 


2010- 


-06- 


■25 


06:19 


15.56 


±0.02 










07:59 


15.57 


±0.03 




2010- 


-07- 


01 


00:13 


15.70 


±0.03 




2011- 


-08- 


-10 


00:41 


15.29 


±0.03 










02:18 


15.41 


±0.03 










03:55 


15.38 


±0.03 




2011- 


-08- 


-13 


00:31 


15.33 


±0.03 










02:06 


15.33 


± 0.02 










23:12 


15.43 


±0.03 


uvwl 


2010- 


-06- 


■25 


06:05 


15.60 


±0.03 










07:48 


15.63 


±0.03 










09:30 


15.61 


± 0.03 




2010- 


-06- 


30 


01:40 


15.59 


±0.02 










03:21 


15.97 


±0.02 










04:49 


15.61 


±0.02 










17:56 


15.71 


± 0.02 










19:32 


15.70 


±0.02 










22:35 


15.77 


±0.03 




2011- 


-08- 


-10 


00:30 


15.36 


±0.02 










02:11 


15.48 


±0.02 










03:50 


15.42 


±0.03 




2011- 


-08- 


■13 


23:01 


15.54 


±0.02 




2011- 


-08- 


-20 


01:00 


15.41 


±0.02 










17:15 


15.38 


±0.04 




2011 


-08- 


■24 


03:10 


15.44 


±0.02 










04:46 


15.45 


±0.02 



Notes. <a) Not corrected for Galactic foreground reddening or extinction 
in the LMC. 
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